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ABSTRACT 

'd\ ' 
o ■ 

O ' We present a flare rate analysis of 50,130 M dwarf light curves in SDSS Stripe 82. 

^ ■ We identified 271 flares using a customized variability index to search ~2.5 million 

^ ■ photometric observations for flux increases in the u- and (7-bands. Every image of a 

^ • flaring observation was examined by eye and with a PSF-matching and image subtraction 

^ ■ tool to guard against false positives. Flaring is found to be strongly correlated with the 

■ appearance of Ha in emission in the quiet spectrum. Of the 99 flare stars that have 

' spectra, we classify 8 as relatively inactive. The flaring fraction is found to increase 

^ ■ strongly in stars with redder colors during quiescence, which can be attributed to the 

• increasing flare visibility and increasing active fraction for redder stars. The flaring 

I ' fraction is strongly correlated with \Z\ distance such that most stars that flare are within 

^ • 300 pc of the Galactic plane. We derive flare -u-band luminosities and find that the most 

^ ■ luminous flares occur on the earlier-type M dwarfs. Our best estimate of the lower limit 
on the flaring rate (averaged over Stripe 82) for flares with Am > 0.7 magnitudes on 

T— I ■ stars with m < 22 is 1.3 flares hour~^ square degree"^ but can vary significantly with 

^ • the line-of-sight. 

cn '. 
O ■ 

^ ■ Subject headings: — stars: flare stars: late-type stars: activity methods: data analysis 

\6 ■ 

o ; 
a^ . 

O ■ 1. Introduction 

> ' 

^ . Flares are a powerful manifestation of stellar magnetic fields, releasing stored energy from 

. surface magnetic loops into the outer atmosphere of the star, and causing strong transient increases 
in the blue and ultraviolet continuum emission as well as in optical, UV and X-ray emission lines. 
There are many unsolved problems in flare physics, ranging from the details of the energy release, to 
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the mechanisms for producing the atmospheric emission, to the understanding of flares on a global 
scale - which stars flare, how often, how much emission do they produce, and on what timescales? 
In the past, the investigation of these global questions has been hampered by the need to obtain 
long time series observations of in dividual stars , which is an exc eedingly time-cons uming exercise 
invol v ing both ground-b ased data (lMoffettlll974l : iLacy et al.lll976l ) and satellite data (lAudard et al 
2OOOI : iGiidel et al.l l2003l ). The advent of large time- resolved sky su r veys c overing a significan t 
fracti on of the sky (e g. Pa n-STARRS, iKaiseil J2004h : LSST, IXvsonl J2002h . Ilvezic et al.l J2008h : 
Gaia. iPerryman et al.l (120011 )) will provide an entirely new opportunity to assemble flare data from 
repeat observations of millions of stars on timescales of minutes to days. The analysis of these new 
data will be quite different from classical flare data which were obtained continuously on single 
objects. 

In this paper we use data from the Sloan Digital Sky Survey (SDSS; York et aD I2OO0I ) to 
perform an initial examination of the possibilities for studying flares in a large sky survey database. 
The SDSS contains repeat observations of one particular part of the sky, known as Stripe 82, along 
the Celestial Equator. This region was targeted because it is in the South Galactic Hemisphere, and 
is visible from Apache Point Observatory during the northern autumn months (Sep-Dec) when the 
main SDSS flelds in the North Galactic Hemisphere are not accessible. The original motivation was 
to obtain numerous photometric scans of the same region of sky so that they could be co-added to 
produce a photometric template 1-2 magnitudes deeper than the normal SDSS data and to quantify 
observational precision and repeatability. It was recognized early on, however, that the individual 
observations formed a valuable time sequence of data that could be mined to identify time-variable 
phenomena. Many additional SDSS-II s cans of this region were also obtained as part of the SDSS-II 



Supernova Surve v (IFrieman et a. 



Data Release 7 (lAbazajian et al 



20081) . The Stripe 82 data products are contained in the SDSS 
2OO9I ) . An investigation aimed at characterizing flare signatures 



in the SDSS Stripe 82 data is therefore timely. 

We have two primary goals in this paper. The flrst is to investigate flare frequencies as a 
function of spectral type and magnetic activity status, where active stars are deflned as those with 
the chromospheric Ha line in emission. We restrict our investigation to low mass stars, of spectral 
type MO and later, and we also concentrate on identifying bona flde flares with high probability 
(rather than accepting lower probability events which may not be flares). Our previous photometric 
and spectroscopic studies of the M dwarfs in SDSS have revealed correlations of the active fraction 
with spectral type (later-type stars are more likely to be active) and with distance from the Galactic 
plane, such that stars clos er to the plane ar e more likely to be active, and statistically are members 
of a younger population (IWest et al.l 12004 |2008| ). The combination of these effects leads to the 
result that later-type stars maintain their activity for a longer time, perhaps for nearly the entire 
~10 Gyr lifespan of the Galactic disk, while early M dwarfs lose their activity in less than 1 Gyr 
( IWest et al.ll2008l ). We wish to determine whether the stars that flare follow the same trends as the 
active stars, as would be expected from the standard model connecting flare emission with magnetic 
reconnection. The SDSS Stripe 82 data provide a unique sample with which to investigate these 
effects, since previous flare data have been obtained only during continuous monitoring of individual. 
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nearby stars. 



The second goal is to provide an initial set of flare frequencies as a function of spectral type 
and Galactic position to compare to our numerical sim ulations of flaring r ates across the Galaxy 
(these simulations are in progress and will be reported in lHilton et al.ll2009l . in p rep). Our previous 
work on the mass and l ur ninosity functions o f the low mass stars in the Galaxy ( jCovey et al.ll2006l ; 
Bochanski et al. 2007bl lal: I Covey et al. 2008b ) together with our knowledge of the active fraction of 
M dwarfs as a function of Galactic height and spectral type, allows us to predict how many stars 
should be flaring at a given Galactic position during a given exposure time. Many of the parameters 
of such a global flaring model can only be roughly estimated at present, and the SDSS Stripe 82 
sample will give us the opportunity to test whether the model is achieving reasonable agreement 
with the data. 

The structure of the paper is as follows. We first describe the SDSS dataset, and our sample 
selection which is designed to choose only those targets with trustworthy observations for analysis 
(§2). We then discuss in considerable detail our efforts to develop a flare variability index that 
can be applied uniformly to the time- resolved photometry to identify candidate flaring epochs (§3). 
The results of applying this flare variability index to the sample are discussed in §4, where we 
investigate the flare rates, correlation with spectral type, activity status, distance from the Galactic 
plane, and line-of-sight through the Galaxy. We also compare our flare luminosities and frequencies 
with previous determinations from nearby stars, and compare the flare rates with simple model 
predictions. In §5, we summarize our results. 



2. Data 



The Sloan Digital Sky Survey (lYork et al.ll2000l : IStoughton et al.ll2002a:|Pier et al.ll2003l:llvezic et al 



Fukugita et al. 



1996 



Gunn et al. 



19981 : 



20041) is a large (~10,000 sq deg), multicolor (uqriz ; 
Hogg et al.ll200ll : ISmith et al.ll2002l : iTucker et al.ll2006l ) photo metric and spectr oscopic survey cen- 
tered on the Northern Galactic cap. The Sloan 2.5m telescope ( Gunn et al. 20061), lo cated at Apache 
Point Observatory, operates in drift-scan mode as the camera (I Gunn et al. I ll998l ) simultaneously 
images the sky in five bands, yielding photometry of ~ 287 million unique objects contained in the 
SDSS Data Release 6 jAdelman-McCarthv et al.ll2008h . 



SDSS Stripe 82 is a narrow band comprising ~300 square degrees along the Celestial Equator 
(—51° < a < +60°, —1.266° < 5 < +1.266°) that was observed more than twenty times under 
photometric conditions between 2001-2005. Many additional (40—60) SDSS-II imaging observations 
of Stripe 82 were obtained du ring 2005-2007; t hese data were not necessarily photometric and 
required additio nal processing (llvezic et al.ll2007l ). The catalog of photometric observations we use 
incorporates the llvezic et al. J2OO7I ) calibrations, and was assembled as a queryable time-domain 
database at the University of Washington in Spring 2008, in advance of the date when the repeat 
Stripe 82 scans were released to the public as part of DR7 (Oct 2008, Abazajian et al. 2008). 
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2.1. Sample Selection 



The following steps outline the process we used to identify a sample of M dwarfs from the 
SDSS Stripe 82 single epoch photometry (DR6), obtain their light curves (from the UW time- 
domain database, though these data are now available in DR7), and restrict the analysis to include 
only reliable, well-determined photometry. Note that in this section and in §3, we perform all 
analysis on the actual measured magnitudes (or fluxes) and their associated errors. In other parts 
of the paper when considering color-spectral type relations, distances, and luminosities, we apply 
interstellar reddening and extinction corrections to the measured magnitudes using the standard 
SDSS correcti o ns co ntained in the photometric database, which are based on the dust maps from 
Schlegel et all JlQQsh . 



1. The SDSS Data Relea se 6 was searched for p oint sources with the co lors of low-mass stars 
in Stripe 82. Based on lWest et al.l (120051 ) and iBochanski et al.l (j2007bl ). the color limits were 
chosen as {r — i) > 0.53 and {i — z) > 0.3 to restrict the sample to M stars. The initial sample 
comprised more than 1.4 million candidate objects. 

2. The photometry in each of the r,i, z filters used to form the colors was required to be brighter 
than the survey limiting magnitudes. This reduced the sample to ~900,000 candidate objects. 



Standard SDSS photometry processing flags (IStoughton et al.ll2002bl ) were used to reject ob- 
jects due to poor DR6 single-measurement photometry. In particular, data were removed if any 
of the following flags were set: SATURATED, NODEBLEND, NOPROFILE, PSF_FLUXJNTERP, 
BAD_COUNTS _ERROR, INTER? _CENTER, DEBLEND_NOPEAK, NOTCHECKED. The 
sample now contained ~700,000 objects from one of the following SDSS Legacy runs: 2583, 
2659, 2662, 2738, 3325, or 3388. 

By far the most significant cut was the next one, where the u photometry was required to 
satisfy a brightness limit u < 22.0. The reason is that flares are most prominent in the u- 
band, and if we cannot measure the quiet (outside of flaring) u magnitude accurately then we 
cannot compare the properties of flaring epochs relative to the quiet star. Because M dwarfs 
typically have (tt — r) > 3, this effectively restricts the sample to stars with r < 19, which 
significantly decreases the volume over which the stars can be observed. Only ~42,000 stars 
passed this cut. 

The UW time-domain database was then employed to search for matches within 2" of the DR6 
position. Observations were kept if they had good astrometry, magnitude errors < 1 mag, and 
good field-wide zeropoint measurements (important for data obtained under non-photometric 
conditions). Only a few hundred stars failed this cut. 

To determine the mean quiet (non-flaring) magnitude in each filter for each star, the measured 
magnitudes were first converted to fluxes (see §3.2 below), and a median flux was calculated 
for each filter. A 3a (where a is the median absolute deviation from the median) cut was 
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then made, and the weighted mean flux (where each measurement is weighted inversely by the 
square of the photometric error) was determined for each filter, subject to the requirement 
that at least five epochs had photometry within 3a available in each filter. Mean fluxes 
were formed for ~41,000 stars. Figure [1] (left) plots the fraction (typically > 60%) of the 
total number of photometric observations that was used to determine the mean flux for each 
star. Figure [T] (right) shows the number of good observations in the u and g filters (that 
is, the observations returned from the UW time-domain database as per step 5 above, but 
not necessarily good according to the SDSS photometry flags). The leftmost asymmetry is 
caused by stars that have u > 21.7 and therefore straddle the detection limit in the w-band. 
The peak near 100 observations is due to the stars that lie in the overlap regions of camera 
columns between the N and S components of Stripe 82, and therefore were observed twice as 
many times. 

7. To ensure that the colors calculated from the light curves are consistent with the initial DR6 
color cut (see step 1 above), the same color cut was applied to the mean quiet magnitudes, 
reducing the sample by a few hundred starsjl] 



8. A color cut requiring that u — g > 1.8 (jSmolcic et al. 2004) was applied to eliminate ~600 



potential unresolved white dwarf-M dwarf binaries, whic h have bluer u — g colo r s. The fre- 



quen cy of M dwarf - M dwarf binaries is likely ~35% fiFischer &: Marcyl Il992l : iReid et al 



20021 ) . Unresolved M dwarf binaries add uncertainty to the analysis, since the color of the 
star may not give an accurate quiescent luminosity and distance. However, it is not possible 
to eliminate unresolved M dwarf - M dwarf binaries using photometry alone. Our final ugriz 
sample consists of 39,633 stars with observations totalling 1,987,253 epochs. 

9. An additional photometric sample was obtained by following the prescription above for stars 
with reliable ug photometry, but riz photometry that is saturated, or faulty in another way. 
Starting with u < 22, u > g, and more stringent magnitude requirements for the u and g 
bands (magnitude errors < 0.1), an additional ~500,000 point sources were obtained from the 
DR6 catalog in Stripe 82. Matching to 2MASS within 2", requiring good i^-band photometry 
(SNR > 5, magnitude error < 0.1) and applying the color cuts g — K > 4 and u — g > 1.8 
leaves ~ 12,000 stars! Matching to the time-domain database within 2", and requiring well- 
defined quiet magnitudes (see step 6, above), gives the ug-only sample of 10,497 stars (524,981 
epochs). 

Although we initially analyzed the two samples separately, the results were very similar and 
showed no systematic bias due to the lack of riz colors for the ■u^'-only sample. Therefore, the 
rest of the analysis is shown for the combined sample. 



^After correcting for interstellar reddening and extinction the limiting r ~ i color was reduced to 0.43 and the 
limiting i — z color to 0.23. 

^After correcting for interstellar reddening and extinction the limiting g — K color was reduced to 3.6. 
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2.2. Spectroscopic Observations and Spectral Types 



About 10% of the stars (4005, comprising 197,867 epochs) in the ugriz photometric sample 
have SDSS spectroscopic observations, determined by matching to the spectroscopic DR6 database. 
None of the ug-onlj sample have SDSS spectra (these stars had poor riz photometry and therefore 
were likely rejected by the spectroscopic targetting algorithm). The stars with SP SS spectra have 
been assigned types using the Hammer spectral typing facility (ICovey et al.l 120071 ) and have been 
previously analyz ed for magnetic activity based on detecting emission i n the chroniospher ic Ha line 
( IWest et al.ll2008l ). Stars characterized as active and weakly active by lWest et al.l (120081 ) comprise 
the active sample (219 stars), with the rest of the spectroscopically observed stars being labeled 
inactive (3786 stars). 

The vast majority of stars in our sample do not have spectra, so we assign an estimated spectral 
type based on the photometric colors. The estimated spectral types serve as a proxy for mass and 
surface temperature and provide a convenient way to subdivide the sample into intuitive "early, 
mid, late" spectral type bins to guide the analysis and compare to M dwarfs which do not have 
SDSS photometry. The colors used to assign types must first be corrected for extinction in order to 
avoid the common problem of assigning a reddened, early-type star at large distance an erroneous 
later spectral type, which would in turn place it at an (incorrect) closer distance. Since Stripe 82 
is located in the Southern Galactic hemisphere, our line of sight to nearly all stars (except the very 
closest) passes through the plane of the Galaxy. We therefore account for the interstellar reddening 
and extinction as described in §2.1. 

The (r — z, i — z) color-co lor diagram w as used to classify t h e ugr iz sample. We fit 2D Gaussian 
distributions using FastMix (IMoord Il999l ) to the IWest et al.l (120081 ) DR5 M dwarf spectroscopic 
sample comprising some 20,000 stars with good riz colors and accurate spectral types based on 
SDSS spectral The resulting best fit parameters (mean and covariance matrices, given in Table [1]) 
define a 2D Gaussian probability function for each spectral type, allowing us to assign estimated 
spectral types based on the photometric colors, as shown in the top panel of Figure [2J 



For the ug-onlj sample, we estimate the spectral types according to the g—K color (ICovey et al. 



2008bl ). Again using the DR5 M dwarf spectroscopic sample, with K magnitudes from 2MASS, a 
ID Gaussian probability function for the extinction-corrected g — K color was determined for each 
M spectral type. The bottom panel of Figure [2] shows the distribution of r — « vs. g — K indicating 
that both colors provide reasonable estimates of the spectral type. 

The dotted line in the top panel of Figure [3] gives the estimated spectral type distribution of the 
photometric (combined ugriz and ug-onlj) sample. The effect of the limiting apparent magnitude 
of the sample [u < 22) is clearly evident through the relatively large numbers of early-type (MO-Ml) 
stars which can be seen to larger distances, and the small numbers of later-type (M4+) stars which 



■^We augmented the original DR5 sample with additional MO stars having bluer colors in order to fit the full color 
range of the extinction-corrected colors in our sample 
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can only be observed nearby. Therefore, the later-type stars come from a much smaller volume than 
the early-type stars. The colored lines in Fig ure Bl show the di stributions of active and inactive stars 
as a function of spectral type. As expected (jWest et al.l 120041 ). most of the early-type M dwarfs are 
inactive, while most of the later ones are active. 

Finally, we investigated the contamination by gi ants in the sanaple. S eparating giant and dwarf 
stars without high resolution spectra is problematic. ICovey et al.l (l2008bl ) examined the J vs. i — J 
color-magnitude diagram of two calibration regions (1 sq deg each) in Stripe 82. We verified that 
the colors of the stars in our sample lie in the region of this color-magnitude diagram with < 2% 
contamination from giants. 



2.3. Followup Spectra from the ARC 3.5m Telescope 

Time domain survey data potentially allows the investigation of flaring rates from both active 
and inactive stars. The former have been well studied for many years, but the latter have received 
very little attention. Because we now have available a sample of "stars that flared" , it is interesting 
to know if they are active or inactive in their quiescent state. Therefore we obtained followup spectra 
with the ARC 3.5m Telescope at Apache Point Observatory, using the Dual Imaging Spectrograph 
(DIS) with the R300 grating giving wavelength coverage ~6000 A-9200 A. On UT080825, we used 
a 2.0" slit giving R ~ 800 and on the dates UT080919, UT081013, UT081107, UT081124, and 
UT081214 we used a 1.5 " slit giving R ~ 1000. As these stars are quite bright in r, exposures of 
3-15 min gave S/N > 10 at Ha. Blue spectra were also obtained; however, the S/N was inadequate 
to obtain information about the activity status of the other H Balmer lines or Ca II H & K. The 
data were reduced using standard IRAF0 procedures, wavelength calibrated using HeNeAr lamp 
spectra, and fiux-calibrated using spectrophotometric standard star observations (though some of 
the data were taken during non-photometric conditions, limiting them to a relative flux calibration). 
We discuss these observations further in the context of the relative probability of flares on active 
and inactive stars in 54.6. 



3. Variability Analysis 

3.1. Properties of Flare Photometry 

The SDSS Stripe 82 photometry for each M dwarf typically consists of ~50 good observations 
(Figured]), irregularly spaced in time with separation of several days, during the fall months, span- 
ning several years. These data allow us to construct a low cadence light curve in each filter. Since 



''IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation 
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flaring timescales range from minutes to liours, individual flares will be observed only once during 
their time-evolution. Therefore, each epoch of each star must be analyzed for a flare signature. An 
example of the ugriz light curves for an active M6 dwarf is shown in Figure HI This star flared 
(with increase^ of 5.5 magnitudes in the w-band, and 3.1 magnitudes in the ^t-band) during the 
first epoch marked with a star and again at much lower levels in the second marked epoch. 

Flare emission can change significantly over the timescale of the SDSS photometric observa- 
tions, which comprise ~54 seconds in each filter as the star passes across the detector, in the filter 
order r-i-u-z-g. For the u and g observations which are primarily of interest here, this means that 
the g observation begins 108 seconds after the u observation, and there is a 54 sec gap between 
them. Figure O illustrates the photometry measurements for a small fiare, with a 0.6 mag increase at 
the peak in the f/-band (similar to u) and a 0.1 m ag increase in the .B-b and (similar to g), obtained 



with the NMSU Im telescope during April, 2008 (iKowalski et al.ll2009l . in prep). Depending on the 



timing of the SDSS observations obtained during such a fiare, our requirement that both the u- 
and (jf-bands show enhancements (as described in the next section) would inefficiently select fiares 
of this size. This is one of the reasons that we limit our fiare detection capability to only those 
events that show 0.7 mag or larger increases in the u-band (see §3.3). 



3.2. Flare Variability Index 

We have developed a Flare Variability Index, ^ug, that is a modified version of an index origi- 
nally designed to rn easure variability in multi-epoch and multi-band photometry (jWelch fc Stetson 
19931 : IStetsonlll996l ). The Welch and Stetson Variability Index is often employed for periodic vari- 
able searches, e.g., for Cepheids and RR Lyrae stars. Our modification produces a variability index 
that is calculated on an epoch-by-epoch basis. The Flare Variability Index, is based on mea- 
suring positive fiux increases in both the u and g filters, since the d efining charact eristic of stellar 
fiares is the dramatic enhancement in the blue continuum emission (jMoffettlll974l ). Note that the 
photometric data in the UW time-resolved database are given in m agnitudes, w h ich w e convert to 



nanomaggies (nMgy), a linear fiux densi ty unit, via Equation 1 in llvezic et a 



2004 



fl2007f). One Mgv 



Blanton et al 



corre sponds to a fiux density of 3631 Jy (lOke fc Gundll983l : iFinkbeiner et al. 
20051 ) . With fiuxes expressed in nanomaggies, we quantify the increases in the u and the g fiuxes 



at a given epoch on a given star by comparing them with the measurement errors: 



\-^u,j ■P'u, quiet) 



ug 



X 



y^aJ Pg,quiet) 



(1) 



where is the Flare Variability Index evaluated at epoch j, F„ , is the fiux in the u filter at epoch 



^Although increases in flux correspond to negative changes in magnitudes, to avoid confusion we refer to flares as 
having positive magnitude changes z\u — \ti^iare '^quiet 

I throughout the paper. 



- 9 - 



j, Fu^quiet is the mean u-finx for the star, and ap^ . is the photometric uncertainty in the u flux at 
epoch j; the same notation is used for the g filter quantities. The quiet flux values Fu^guiet and 
Fg^quiet are taken to be the weighted means of the well- measured epochs (see §2.1). 

Flares will have positive flux deviations in both u and 

Fu,j Fuquigl > 0, Fg j Fgqy^^gf > 0, (2) 

generating a positive value of Note that eclipses, which have negative flux deviations in both 
filters, also generate a positive value of while the subset of random noise characterized by 
positive flux deviation in one filter, negative in the other, generates a negative value of To 
insure that eclipses are not included as candidate flare events, we require that the w-band flux 
enhancement be positive. 



3.3. Flare Candidate Selection 

The quiet flux of M dwarfs in the n-band is very faint, in most cases close to our apparent 
magnitude limit ofu = 22. Thus we expect that photometric variations caused by random sampling 
of noisy data (e.g. the significant noise apparent in the u and g light curves in Figure HI) will in 
some cases produce positive values of the Flar e Variability I ndex, and hence masquerade as flares. 



The False Discovery Rate (FDR) method of iMiller et al.l (120011 ) allows us to further limit the 
number of positive $ng events that are candidate flares. The FDR analysis involves comparing 
a source distribution with a null distribution, and then using the null distribution to assign a 
probability that a given value of is drawn from the null. Figure M compares the candidate flare 
(source) distribution {^ug > 0, F„j — F^^quiet > 0, Fgj — Fg^quiet > 0), with the random noise (null) 
distribution {^^g < 0), for the same number of epochs in each distribution (~770,000). It is clear 
that a significant number of the candidate flares, especially at low values of ^ug, can be attributed 
to random noise. False Discovery Rate analysis provides a way to select a threshold value that 



limits the percentage of false positives. Employing the IDL routine in Appendix B of iMiller et al. 



( l200ll ). we determined that setting the FDR parameter, a, to 10% corresponds to ~ 100 (the 
vertical line in Figure [6]). Therefore, no more than 10% of the epochs with ^ug > 100 are false 
positives caused by random sampling of noisy data. This returned a reasonable number (673) of 
candidate flare epochs for individual examination. 

As shown in Figure O and described briefly in §3.1, flares with small flux enhancements typically 
occur over short timescales and may not be detected in both the u and g filters during the SDSS 
imaging cadence. An additional complication is that the SDSS u-hand imaging suffers from a 
photometric defect known as the "red leak" caused by the u filter which allows flux with wavelength 
longer than 7100A to be partially transmitted. Due to the complexity and time dependence of this 
instrumental artifact, it is not corrected for in the SDSS photometryO We investigated the red leak 



see http:/ /www. sdss.org/dr7/products/catalogs/index.html for a summary description of the red leak 
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by analyzing SDSS images taken at a range of airmasses. The differential refraction between the 
blue and red light separates the red leak photons into a faint image that is offset from the (blue, u- 
band) stellar image at the pa rallactic angle and with the correct image displacement, ranging from 



1.15"to 2.4^'f lFilippenkdll982l ). As a result, there may be various degrees of PSF deblending from 
observation to observation depending on the seeing, airmass, and stellar and flare colors. We show 
examples of the red leak effect in Figure [71 In the two worst cases, we estimate that the amount 
of the red leak would have added ~0.3 and 0.5 mags to the w-band measurements had they been 
co-aligned and therefore not deblended. Thus, we do not trust any measured w-band magnitude 
enhancements of 0.5 magnitudes or less. Since, as described in the Introduction, we only want to 
analyze bona fide flares, we adopt a conservative limiting magnitude of Au = \ufiare — Uquiet\ > 0.7 
magnitudes, which corresponds to almost a factor of two increase in the flux, as the lower limit of 
our flare detection capability. Requiring a correlated increase in g is an additional safeguard against 
erroneously selecting "flares" due to the red leak. 

Flares typically show blue conti n uum emission with the approximate spectral shape of a 



10,000K blackbody (IHawley fc Fished (jl992l ). see also Figure [5]). This characteristic implies that 



the flux enhancement in u should be greater than that in g, 



F ■ F 

> = eg (3) 



u, quiet g, quiet 

All but two of the candidate flares exhibit this characteristic. The time evolution of the flares 
combined with the image sequencing adds ambiguity to the precise spectral shape of the flares, 
and so we do not exclude candidate flares that fail to exhibit this characteristic. In rare cases, 
observations during the rapid rise phase of the flare, if the «-band was measured just as the flare 
started, and the gf-band slightly later, at the peak of the flare, might show > Cu- 

The number of candidate epochs is now small enough that each can be queried in the time- 
domain database to return photometry flags, as was originally done for the SDSS DR6 Stripe 82 
photometry (see §2.1). The following flags cuts were applied in the u and g bands to eliminate obser- 
vations with poorly determined photometry: SATURATED, NODEBLEND, NOPROFILE, PSF_ 
FLUXJNTERP, BAD_COUNTS_ERROR, INTERP_CENTER, DEBLEND_NOPEAK, NOTCHECKED, 
MAYBE_CR and SATURATED = 0. Note that most of the rejected epochs resulted from having 
NODEBLEND set, indicating that a star is seen as two fainter objects in good seeing, but a brighter 
object (masquerading as a flare) in poor seeing, when the two objects cannot be deblended. 

The candidate flare sample has now been pared as much as possible by statistical analysis, 
noise requirements, and flare physics. Further analysis requires visual inspection of the data for 
each epoch. An image was retrieved and examined to look for artifacts such as diffraction spikes. 
The light curves were also examined to insure that the data are consistent with those seen from 
flares (e.g. Figure H]). Finally, to confirm the increases in flux indicated by the automatic SDSS 
PSF photometry, we identified a non-flaring epoch for each star that had a candidate flare epoch, 
and ran a PSF-matching and image subtracti on algorithm develo ped for the SDSS-II Supernova 



Survey and the Deep Lens Transient Survey (iBecker et al.ll2004j ). During a true flare, the field 
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stars will disappear in the subtracted images, and the flare star will show residual PSF structure in 
both the u and g images. For cases with unusual PSF structure or data taken with poor seeing, we 
calculated the aperture photometry for the flare star and compared to several other stars in the field, 
as the aperture photometry should be independent of the PSF shape. Only 9 flare candidates were 
eliminated through this visual inspection, image subtraction, and aperture photometry process. 

Figure [H] shows the selected flare candidate region in a plot of w-band flux enhancement versus 
Flare Variability Index. Black points show all observations with $„g > 0, Fuj—Fu^q > 0, Fgj—Fg^g > 
0. The vertical line corresponds to = 100 and the horizontal line corresponds to e„ = 1.91 (Am 
mag = 0.7). The 271 flare epochs are shown as circles, colored by most probable spectral type. To 
the left of our threshold, observations may have a large flux enhancement; however they are 
not selected due to either a large photometric error in m or a negligible increase in g. We inspected 
several candidate epochs that fell just short of making our cuts, and found that the return in actual 
flare observations was very small compared to the increase in the number of candidate epochs 
resulting from weakening our criteria. 

3.4. Summary of Flare Sample 

Table [2] summarizes the steps we took to arrive at our final flare sample, which consists of 271 
epochs (225 from the ugriz sample and 46 from the ug-only sample) on 236 stars. We verified that 
the fraction of observations that flare for all the runs is fairly uniform, consistent with observing 
random stellar events and not the result of photometric artifacts on one or a few nights. There 
are 38 flares (38 stars) on early MO-Ml dwarfs out of 1,783,142 epochs (35,024 stars) observed; 83 
flares (79 stars) on mid M2-M3 dwarfs out of 629,170 epochs (12,986 stars) observed, and 150 flares 
(119 stars) on later M4-M6 dwarfs out of 99,922 epochs (2,120 stars) observed. Clearly the flaring 
fraction is much higher on the later-type stars, both because they are more likely to be magnetically 
active, and because flares at the level of twice the quiescent brightness are much lower energy on 
fainter stars, and therefore likely to be observed more often. These results are discussed further in 
§4. II below. Note that several of the redder stars flared multiple times, with the highest occurrence 
being 5 flares (and possibly several more small flares that failed to make our cuts) on an active M5 
star (SDSS J025951. 71+004619.1). 

4. Results & Discussion 

4.1. The Fraction of Stars that Flare 

The Stripe 82 imaging data provided observations of 2, 512, 234 epochs with reliable photometry 
in both the u and g bands, from which we confidently identified 271 flares. This gives an overall 
flaring fraction of 0.0108 ± 0.0007%. In other words, about 1 of every 10, 000 observations in 



-VI- 



Stripe 82 of a star with the colors of an M dwarf and m < 22 resulted in a flare detection with 
Aii > 0.7 mags. Several factors enter into this determination of the flaring fraction: the stellar 
density distribution being sampled to an apparent magnitude hmit of < 22, resulting in a large 
number of early type M dwarfs in the sample; the fact that active stars are more likely to flare, and 
are more likely to be of later spectral type and therefore are being sampled over a smaller volume; 
and the flare detection limit of Am > 0.7 magnitudes which results in lower energy flares being 
detectable only on intrinsically fainter stars. We discuss each of these effects separately. 

As expected, the stars that were known to be active, based on SDSS spectral observations 
showing Ha in emission, flared much more frequently. Out of 10,396 observations of known active 
stars, we found 29 flares, which gives a flaring fraction of 0.28 ± 0.05%, ~30 times higher than the 
flaring fraction for the whole sample. There were no known (from SDSS spectra) inactive stars that 
flared (but see §4.2 below). As far as we are aware, these are the first determinations of the flaring 
fraction of M dwarfs (and active M dwarfs) as a Galactic population. 

Redder stars (later spectral types) also were found to flare more frequently than earlier type 
stars, as expected because they are more likely to be active and because lower energy flares are 
detectable due to their lower surface temperatures and fainter quiescent fluxes. The top left panel of 
Figure [9] shows the flaring fraction separated by spectral type group. The fractions differ by about 
an order of magnitude for each group, ranging from 2 x 10^^ for the MO-Ml stars to 1.3 x 10~^ 
for the M2-M3 stars to 1.50 x 10^^ for the M4-M6 stars. It is notable that the M4-M6 stars had a 
factor of 18 fewer observations (100,000 compared to 1.8 million) yet had nearly 4 times as many 
flares. Clearly the later-type stars form the bulk of the flare sample. 

Figure [3] shows that the total sample is sparsely populated by stars with spectral types M5 and 
later since they must be nearby to pass the -u < 22 threshold. In fact, ~96% of the observations 
in the M4-M6 bin are of M4 stars. If we examine the M4 and M5-M6 bins separately, the flaring 
fractions become 1.2 x 10~^ ± 1 x 10"^ for M4's and 8.5 x 10"^ ± 1.5 x 10~^ for M5-M6's, continuing 
the order-of-magnitude increase in flaring fraction towards later-types. However, the small number 
of stars (82), epochs (3656), and flares (31) on M5-M6 stars in the sample makes this value relatively 
uncertain, so we don't include it in Figure O 

If we remove the dependence on activity by plotting the flaring fraction of only the known 
active stars for each spectral type group (bottom left panel of Figure [9]), the fractions are signifi- 
cantly higher but the trend with spectral type remains. The detection limit. Am > 0.7, may be a 
contributing factor to the high flare fraction on the later spectral type stars. We can investigate 
this effect by considering only large flares, where a large flare is defined as having a -u-band lumi- 
nosity greater than 1% of the bolometric luminosity (see §4.4 for details on the calculation of flare 
luminosities). These higher luminosity flares will be detected at all spectral types in our sample. 
Although there are fewer flares that are large (73), they provide a comparison of the flare fractions 
between early-, mid-, and late-types without a threshold bias. Figure [HI (bottom right panel) shows 
the fraction of flaring epochs with large flares. The flaring fraction for the later-type M dwarfs is 
still an order of magnitude larger than for the early-types, indicating that the detection threshold 



is not the main factor contributing to this result. 



4.2. The Flaring of (Relatively) Inactive Stars 

In the SDSS spectroscopic sample, every star that flared showed Ha emission in its quiescent 
spectrum. To further investigate the activity status of the stars that flared but had no SDSS spectra, 
we obtained foUowup spectra of 76 additional stars using the DIS spectrograph on the ARC 3.5m 
Telescope (see §2.3). Thes e spectra were spe ctral typed and had Ha emission measured using the 



Hammer analysis software (jCovey et al.ll2007l ). The results are given in Table O Seven stars show no 
evidence of Ha emission. One star, SDSS J022112. 48+003225. 7, has an ambiguous signature near 
6563 A, but this feature may be a blend of Ha and a TiO band known to exist in inactive M4 stars 
at this wavelength. Without higher resolution spectra, we group this star among the inactive stars. 
Another star, SDSS J204452. 77+010546.9, exhibited variable Ha emission between no activity and 
weak activity on two different epochs (UT081107 & UT080919). Perhaps we happened to catch the 
UT080919 spectrum at the tail end of a flare; however, an analysis of the low-signal B400 spectra 
(3400 A— 5500 A) showed no significant change in color between the two epochs, indicating that 
the observation would have likel y not met our fla ring criteria. We classify this object as variable, a 



known trait of active M dwarfs (IGizis et al.ll2002l ). 



Of the eight likely inactive stars, four were observed on two separate nights to confirm the lack 
of emission. The spectral types for the eight stars show no obvious trend toward early- or late- 
type (Hammer type (photometric type): K7 (MO), MO (MO), MO (MO), Ml (MO), Ml (MO), 
M3 (M3), M4 (M4), M4 (M3)) but the numbers are too small to draw any conclusion about the 
spectral type distribution of inactive stars that flare. Though formally classified as inactive, these 
stars may b e part of a group of weakly a ctive M dwarfs which show Ca II H and K emission but Ha 



absorption (jWalkowicz fc Hawleyl 120091 ). Certainly the fact that they flared is a strong indication 



that they possess surface magnetic fields. 

From our combined spectroscopic sample (SDSS + DIS) of (23+76 = 99) stars, with 8 being 
classified as inactive due to lack of Ha emission, we derive a preliminary fraction of flaring stars 
that are inactive to be ~8% (the binomial errors on this percentage are 1.9%, +3.6%). In other 
words, out of every 12 stars that flare in the SDSS Stripe 82 region, one of them might be inactive. 
This is an important value to determine for modeling the rate of flares in the Galaxy, since as noted 
in §2.2, most M dwarfs are inactive, especially at early-mid spectral types. 



4.3. The Flaring Fraction vs. Distance from the Galactic Plane 

As a result of dynamical encounters, stars can move farther from the Galactic plane as they age. 
The mean stellar age thus tends to increase as a function of vertical distance \Z\ from the plane. 
Therefore, other stellar properties that change with the lifetime of the star should also change 
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with \Z\ distance. The fraction o f active stars in the SDSS DR5 sample was found to decrease 
with \Z\ for each M spectral type (IWest et al.ll2008l ). indicating that magr ietic activity depends on 



age. Many o t her studies have linked chromospheric activity to stellar age (IWi 



Wielenlll977l : iGiampapa &: Liebertlll986l : ISoderblom et al.lll99ll : iHawley et al 
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19961. 
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Figure [TO] shows the dependence of the flaring fraction on vertical distance (age) for our three 
spectral type groups. Computing the fraction of epochs that flare effectively normalizes by the 
local stellar density and accounts for varying incompleteness i n different distance bin s. Distances 
were estimated using an r — 2; photometric parallax relation (IBochanski et al.ll2009l. in prep) for 
stars in the ugriz sample and a g — K photometric parallax relation (jCovey et al.ll2008al ) for stars 
in the it^r-only sample. As discussed previously, to obtain accurate estimates for the distances we 
applied interstellar extinction corrections to the r — z and g — K colors before using the photometric 



parallax relations. The S un was taken to be 15 pc above the Galactic plane (ICohenlll995l : iNg et al 



19971 : iBinney et al.lll997l ). The fraction of epochs where flares were detected decreases sharply with 
\Z\ distance for all three spectral type groups. This is expected because of the age dependence 
of (quiescent) magnetic activity, and the strong correlation between active stars and stars that 
flare. As mentioned in §2.2, the volume over which the stars are measured differs markedly for the 
different groups. The later-type group is sampled only at close distances, while the mid- and early- 
type groups are sampled progressively farther from the plane. The cumulative distributions for the 
flaring epochs (asterisks), active epochs (triangles), and all epochs (open squares) are given in the 
bottom panel of Figure (TU] and show the combined effects of age, spectral type and stellar density 
distribution in an apparent-magnitude limited sample. More than half the stars in the sample are at 
a distance greater than 300 pc from the plane, but less than 10% o f the flares occur on stars that far 
away. Since the thin disk outnumbers the thick disk until ~1 kpc (jjuric et al. 2008), it is likely that 
these distant stars that flare still belong to the thin disk population. Although 70% of the active 
epochs are within 300 pc, a third of the active stars are early-type stars which do not contribute 
significantly to the flaring population. For these reasons and because of higher quality photometry 
for stars with u < 21, the nearby late-type, active stars clearly make the largest contribution to the 
flaring fraction in this sample. 

The flaring f ractio n shows a larger decrease with increasing \Z\ distance than the active fraction 
from lWest et al.l (120081 ) . especially for the M4-M6 bin. We find that an additional factor contributing 
to the trends seen in Figure [10] is a bias towards higher l uminosity flares on more distant stars. M 
dwarfs are known to exhibit a large spread m u — g color (jCovey et al.ll2008a] ). For a given spectral 
type, we find a large spread in w-band quiescent luminosity (see §4.4 for details on the calculation of 
flare luminosities). For example, the M4's in this sample have an average ti-band luminosity during 
quiescence of ~4xl0^'^ erg s~^ with a spread of ~2xl0^^ erg s~^. For distances > 200 pc, most M4's 
become too faint in the n-band to be included in our sample. At these large distances, only the 
M4's which are intrinsically bright are included (i.e., there is a Malmquist bias due to our apparent 
magnitude selected sample). Therefore, a flare with Au = 0.7 magnitudes on an intrinsically 
bright M4 has a larger w-band luminosity than a flare with the same apparent enhancement on 
an intrinsically fainter M4. More luminous flares (or more luminous stages of flares) occur less 
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frequently (ILacy et al.lll976l ). which generates a bias in our sample towards smaller flare fractions 
at greater distances. A useful future study would be to examine flare rates for a flare-luminosity 
limited sample. 



We also note that IWelsh et al.l ( 20071 ) found that most of the M dwarf flares detected in NUV 
Galex observations (1750 A— 2750 A) occur on stars within 300 pc of the plane; the flare flux in 
the NUV- and u-bands is thought to come from the same (or similar) emission mechanism, so it 
is expected that flare detections in these wavelength regimes should trace one another through the 
Galaxy. 



4.4. Flare Magnitudes and Luminosities 



Figure [TTl top, shows the distribution of flare magnitude (An) for the three spectral type 
groups. The later-type stars generally show larger magnitude increases, due to their lower quiescent 
luminosity. In other words, it is much easier to see an energetically small flare on a star that is 
intrinsically fainter. Since the flare detection limit was chosen to be Au > 0.7 mags on all stars, 
this corresponds to a lower luminosity fla re on the intrinsically fainter star. T his relationship is 
consistent with previous flare studies (e.g.. lMoffettlll974j : iRockenfeller et al.ll2006l ). Figure [TT] shows 
empirically which stars are likely be observed with large magnitude increases, which is the metric 
of interest for time domain surveys searching for transient brightness variations. 

To obtain flare luminosities, the quiescent stellar -u-band flux was subtracted from each flare 
observation, and the fluxes (which are per Hz) were converted to luminosities using the distance 
to each star (see §4.3) and the nominal u bandwidth (FWHM 1.33x10^"^ Hz, or roug hly 600 A). 
Again, we use extincted-corrected magnitudes for the luminosity calculations. The flare luminosity 
distributions by spectral type are illustrated in the bottom panel of Figure [11] We calculate the 
errors on the luminosities to be ~50%, which is dominated by the error in the photometric distance 
determination and includes the uncertainty in the amount of red leak flux (up to 30%) contained 
in the measured quiet u magnitude. Systematic errors in the luminosity calculation are introduced 
since we assume full line-of-sight corrections for interstellar extinction. As mentioned in §2.2, 
corrections effectively increase the distances, thereby increasing the luminosities by ~10-70%. 

In Figure [TT], the low luminosity cutoff of each distribution represents the detection limit 
of A-u = 0.7 mags, which corresponds to nearly an order of magnitude difference in luminosity 
for each group (log L ~ 28.0 for M4-6, ~ 28.7 for M2-3, ~ 29.5 for MO-1). Although the flare 
magnitude enhancements on the earlier-type stars are smaller (as shown in the top panel), their 
derived luminosities are genera l ly lar ger t han the fla r es oii the later-type stars. This is consistent 
with the flndings of iLacy et al.l (119761 ) and jPettersenl (11988! ) . that the average rate of energy loss in 
the f7-band during flares increases as a function of stellar bolometric luminosity. For example, if we 
compare the two largest -u-magnitude enhancements in the SDSS flares, the Au = 4.7 flare on a Ml 
star gives a luminosity of 1.3 x 10^^ erg s~^ while the Au = 5.5 flare on a M6 star gives a luminosity 
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of only 1.7 X 10^^ erg s^^. However, we emphasize that our derived flare luminosity is sampled at 
a random time during the flare and there is no way to determine whether the emission occurred 
during the rise, peak or decay phase so it is not strictly comparable to the average energy. It is 
interesting that we do not observe very high luminosity (> 10^° erg s~^) flares on the very red stars. 
This is likely due to two factors: they are quite rare on the relatively low (quiescent) luminosity 
stars; and the number of epochs that SDSS observed on the later-type stars is much smaller than 
on the earlier type stars since the latter comprise a much larger fractioi i of the sample. Whe n we 
compare the flare luminosity to the bolometric luminosity of the stars ( Reid fc Hawleyl [2OO5I ) . we 



find a few flares that emit as much as ~1-10% of the bolometric luminosity in the w-band, on stars of 



both early and late spectral types. iPetterseru (119881 ) had determined a limit of 10 for the average 



Ljj / Lf)oi] however, our values may represent instantaneous measurements of more luminous stages 
during flares. 

With a common flare detection threshold {Au = 0.7 magnitudes) for all spectral types, the late- 
type stars emit less luminous, yet many more, flares. However, we find that the total flare energy 
emitted on early-, mid-, and late-type stars is approximately equal for equal observing times. For 
example, the peak luminosities of Figure [TT] occur at log L ~ 30.5, 29.3, and 28.7 erg s~^, a ratio of 
~ 60 : 4 : 1 for early-, mid-, and late-types, respectively. When we sum the luminosities, multiply 
by the filter integration time, and scale to the same total number of observations (§3.4), the total 
flare energies occur in the ratio ~ 0.7 : 0.3 : 1. Apparently, the larger flare rate of the late-types 
approximately makes up for the lower energy output per flare, suggesting that early- and late-type 
stars, summed over a range of quiescent activity levels and stellar ages, may possess similar amounts 
of total stored magnetic energy that is eventually released as flares, but that this flare energy is 
released at different rates on stars of different masses. Further study of a sample including smaller 
flares is needed to confirm this preliminary result. 



4.5. Flaring Frequency 

A class i c stud y of flare frequency for a handful of nearby, very active M dwarfs was carried out by 



Lacv et al.l (Il976f) . They used a monumental collection of continuous photometric monitoring data 
(jMoffettlll974l ). which provided information throughout the time evolution of each flare. An average 
luminosity for each flare may be calculated by dividing the reported total t/-band energy of the flare 
by its duration. In Fig. [T21 we compare the luminosities of the SDSS flares (sampled at a random 
time during the flare) to the average flare luminosities of the Moffett flares for the dM4.5e star YZ 
CMi. In order to estimate the frequency with which a flare of a particular luminosity (or greater) 
was seen, we used the total -u-band observing time for active stars with SDSS spectra and spectral 
types of M4 or M5 (43.335 hours), and obtained a cumulative distribution of flare luminosities. As 
shown in Figure [T^l our derived flaring frequency from the SDSS data is remarkably consistent with 
the flaring frequency of YZ CMi, following a power law form with more luminous flares occurring 
less frequently. 
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4.6. Spatial Flare Rate in Stripe 82 

A number of interest to characterize transient events in large sky surveys is the spatial flare 
rate, the number of flares per hour per square degree. Of course this rate will vary depending on 
which area of the sky is being observed (e.g. with Galactic latitude). For the Stripe 82 data, which 
span / = (45° to 191°); b = (-23° to -64°) the rate can be empirically calculated by taking the 
number of observed flares (271), and dividing by the spatial extent of the Stripe 82 time domain 
database (279.7 square degrees) and the average observing time (50 epochs x 54 sec per epoch = 
2700 sec). This gives a rate of 1.3 flares hr~^ sq deg~^ with a u mag increase of at least 0.1 mags. 

The derived spatial flare rate is a lower limit to the intrinsic flare rate because we have an 
inefficiency, e, for detecting flares. Major contributions to e are varying contrast levels for the 
different spectral types, the fact that not all flares produce an increase in the u and g bands (e.g., 
flares that begin after the u but before the g observation), our limit of 10% false positives by using 
a $„g threshold of 100, and our strict u < 22 magnitude cutoff. As noted explicitly, we intentionally 
adopted a conservative Au = 0.7 magnitude limit for flare detection, so our spatial flare rate is only 
representative of relatively bright flares. 

We calculated the number of flares divided by the average observing time per star (~2700 sec) 
for every square degree (in RA and DEC) in Stripe 82 to investigate how the number of flares hr~^ 
sq deg~^ varies with the line-of-sight through the Galaxy. We find a maximum flare rate of 8 flares 
hr~^ sq deg~^ occurring closer to the plane (/ ~ 56°, b ~ —38°) where the density of stars is higher 
(~240 stars sq deg~^ with w < 22 compared to the average for Stripe 82, ~180 stars sq deg~^ with 
u < 22) and where a higher fraction of stars are active (since they are closer to the plane and 
therefore likely to be younger). 

To illustrate how the spatial flare rate varies with Galactic latitude and longitude, we calculated 
the average flare rate for 3 latitude regions and 3 longitude regionsJil Our results are presented in 
Table HI This raw spatial flare rate decreases with increasing Galactic latitude, as expected due 
to the decreasing stellar density farther from the plane. To account for stellar density, we divided 
by the number of total epochs in each region and present the flaring fraction as a function of 
Galactic latitude and longitude in Figure [131 There is no obvious trend with longitude; however, 
the flaring fraction generally decreases with increasing latitude, probably because lines-of-sight at 
higher latitude include a smaller fraction of stars at small \Z\ distance (see §4.3). 

Finally, the empirical flare rate for all M dwarfs would also include the flares occurring on the 
M dwarf - white dwarf binary systems that we excluded in step 8 of §2.1. We have searched those 
light curves and find 6 events that would satisfy our flare criteria. This has a negligible effect on 



^When binning Stripe 82 into {l,b), we divided the flare rate in each 1 deg x 1 deg square by cos(6) in order to 
correct for the changing area due to converging hues of longitude. To compensate for incomplete coverage at the 
edges of Stripe 82, we also divide by the fraction of each 1 deg x 1 deg bin that was scanned by Stripe 82, and we 
consider only bins that have >50% coverage. 
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the overall flare rate. 



4-6.1. Comparison to Simple Model Prediction 

For comparison, we can also estimate the flare rate from our knowledge of the spatial density 
of M dwarfs in Stripe 82 (Figure 3), the fraction of active stars at each spectral type (Figure 3), the 
relative flare rates for active an d inactive stars (^ 4.1, 4.2), and the flare frequency distributions for 



stars of different spectral types (jLacy et al.lll976l ). This simple model calculation gives a prediction 



of the spatial flare rate we should observe in Stripe 82. 

We consider the same spectral type ranges, and the values measured or derived previously to 
form the general equation: 

# flares = S X {f active X fajlare + f inactive X fijlare) X 1/ X t (4) 

where s is the number of stars, f active is the fraction of stars that are active, fajiare is the fraction 
of active stars that flare, finactive is the fraction of stars that are inactive, fijiare is the fraction 
of inactive stars that flare, u is the number of flares hr~^ above a limiting luminosity, and t is 
the observing time per star. A very uncertain parameter is fijiare, and we calculate the predicted 
number of flares adopting both 1% and 10%]^ We list the measured and derived values for Equation 
4 in Table E 

In §3.4 we found that the number of conflrmed (observed) flares was 38 for MO-Ml (compare 
to 20 - 66 predicted); 83 for M2-M3 (compare to 68 - 150 predicted); and 150 for M4-M6 (193 - 
233 predicted). In all cases the predicted values are reasonably close to the observed number of 
flares. Given the crude estimates in our predictions, this agreement is encouraging, indicating that 
our method for flnding flares in the low-cadence SDSS photometric data is apparently reliable. 

There are also several places where the predictions can be improved. Two very poorly de- 
termined parameters are the fraction of inactive stars that flare and the flaring frequency, Ui, of 
inactive stars that flare (here taken to be the same as the flaring frequency, v, of active stars of the 
same spectral type). We have an ongoing pr ogram to investigat e the flaring rates of inactive stars 



through continuous monitoring observations (IHilton et al.ll2009l . in prep). 



Furthermore, our estimates for ly are taken from YZ CMi, a star w ith a well-determined but 



likely higher-than-average flare frequency distribution (iLacy et al.l 119761 ). Therefore, it is not sur- 
prising that our predicted number of flares for this bin exceeds the observed number. Even if we set 
I'i = flares hr~^, 189 flares are predicted. If YZ CMi accurately represents the flaring frequency 
for active M4-M6 stars, we estimate our inefficiency, e, for detecting flares above A-u > 0.7 mags on 



^In §4.2 we found that ^^8% of the flare stars are inactive. Although this is not strictly the same as fijiare, which 
is the fraction of inactive stars that flare, the 8% value made the adopted range of 1 - 10% seem reasonable. 
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M4-M6 stars to be > 20%, assuming all the other parameters are correct. Alternatively, if instead 
we use a slightly lower flaring frequency for all M4-M6 stars, v = 0.3 flares hr~^, we predict 145 — 175 
flares, which is consistent with the observed value. 

We also note that the active fraction of M4 stars found here (30%) is considerably lower than 
the active fraction found from the full DR5 sample (50%, West et al. 2008). DR5 is dominated 
by stars in the Northern Galactic cap, with most observations obtained looking out of the Galactic 
disk. In contrast. Stripe 82 is observed through the disk, looking toward the Southern Galactic cap. 
The difference in active fraction between the North and South is intriguing but unexplained. We 
are pursuing an additional study of the change in active fraction with galactic latitude. 

Several regions in Stripe 82 have moderately high extinction (> 0.5 mags in the g-hand), and 
we examined the possibility that the spatial flare rate is correlated with star-forming regions and 



therefore contaminated by accretion events. We compared the composite Ha map from iFinkbeiner 



(120031 ) to the spatial distribution of stars that flare and note an absence of any conspicuous clustering 
of flare stars near Ha regions. The star-forming regions MBM 18/MBM 19, found at RA > 55°, 
coincide wi th the areas of greatest extinction in Stripe 82 and are suspected sites of low-mass star 



formation (IMcGehed l2008l ) . None of the 11 flare stars with RA > 55° have spectra to examine for 
accretion features. Of the stars that do have spectra, none have especially high equivalent widths 
characteristic of accretors. 



5. Conclusions 

We used the multi-epoch ugriz photometry of the SDSS Stripe 82 to characterize the M dwarf 
flaring properties as a function of spectral type, level of magnetic activity, distance from the Galactic 
plane, and line-of-sight through the Galaxy. In contrast to studying flare rates on a few individual 
nearby stars via continuous monitoring, we quantify the flaring properties of populations of M 
dwarfs using low-cadence observations. We search for a flare signature on an epoch-by-epoch basis 
over millions of observations and employ the FDR algorithm to limit the number of false positives. 
We find 271 flares among ~2.5 million observations. The fraction of epochs that flare increases by 
nearly an order of magnitude for each subsequently redder M dwarf color bin. We attribute this 
trend to there being better 'flare visibility' and an increasing fraction of magnetically active stars 
in the redder bins. 

We find a strong correlation between stars that have flares and stars that have Ha in emission 
during quiescence. For the sample as a whole, 1 out of 10,000 observations shows a flare, but ~30 
out of 10,000 epochs observed on active stars show flares. This finding is consistent with a simple 
physical picture connecting magnetic activity to flaring: surface magnetic loops (the source of the 
persistent Ha emission line) reconnect and accelerate charged particles, which then deliver energy 
to the lower atmosphere and produce the optical flare emission. However, the detailed physics of 
how such large energies in the w-band are generated remains a mystery. Unfortunately, the time 
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difference between the u and g filters adds ambiguity to the interpretation of the flare emission 
spectrum and prevents us from constraining the broadband spectral shape of flares with these data. 

Our total sample is dominated by inactive (no Ha emission) M dwarfs, and we find only eight 
flaring stars that are classified as inactive. Inactive stars may have a much smaller fiUing factor 
of magnetic loops such that they do not show quiescent activity but still have a (small) chance of 
producing flares. Continuous monitoring of inactive M dwarfs is needed to more accurately quantify 
the flaring rates and properties of these stars. 

The amplitudes and luminosities of the SDSS flares are consistent with previous observations. 
Most flares (80%) have amplitudes less than 2 magnitudes in it, which is not surprising since random 
imaging is most likely to catch a flare during its more extended, smaller-amplitude decay phase. 
Both early- and late-type stars exhibit w-band magnitude enhancements as high as Atx ~ 5; however, 
the actual luminosity for a given magnitude enhancement is dependent on the intrinsic shape of 
the underlying stellar spectrum. Very large luminosity flares are not seen on the reddest of stars, 
which may be a combination of a real physical effect and a result simply from there not being as 
many observations of very red stars in the sample. The most luminous flares in the sample emitted 
~10% of the stellar bolometric luminosity at the time they were sampled. 

The flaring fraction is found to strongly decrease with increasing \Z\ distance for all spectral 
types, as expected, due to a smaller active fraction and a Malmquist bias at greater distances. Our 
best estimate of the lower limit on the flaring rate (averaged over Stripe 82) for flares with Au 
> 0.7 on stars with it < 22 is 1.3 flares hour~^ square degree"^, and the larger spatial flare rates 
and fractions are typically found along lines-of-sight at lower Galactic latitude in Stripe 82. Future 
time-domain surveys, such as LSST, will be able to probe to a much larger distance from the plane; 
however, we predict that the spatial flare rate for flares of Au > 0.7 mags on stars from MO - 
M4 will not increase appreciably for intermediate to high Galactic latitudes due to our observed 
\Z\ distance-flaring relation. For lincs-of-sight through the Galactic plane, there is likely to be a 
significant increase in the observed flare rate. Since the ti-band will be observed only infrequently 
with LSST, the g- and r-bands will be the most useful to identify flaring epochs. We find that 
~95% of the flares have g magnitude enhancements of < 1 magnitude. However, due to the higher 
precision (~5 mmag at the bright end; Ivezic et al. 2008) of LSST, a larger number of small 
amplitude flares will be observed. For these reasons, flares will likely form a larger population of 
transients in LSST compared to SDSS, which was limited by the photometric uncertainties in the 
M-band. Since the SDSS flares were selected from the u- and (ji-bands, the properties of these flares 
in the other bands could be applied to candidate flares in future surveys which do not use the 
It-band. 
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Fig. 1. — Left - The fraction of epochs used for the quiescent magnitude calculation shows that 
a significant number of observations were excluded due to being > 3 median absolute deviations 
away from the median. Right - The number of good epochs per star in both the u and g filters. 
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Fig. 2. — Top - The extinction-corrected r — z vs. i — z color-color diagram for the sample of 39, 633 
M dwarfs with well-defined light curves. The spectral types were probabilistically assigned based 
on our presented spectral type-color relations (Tabled]). Bottom - The extinction-corrected r — i 
vs. g — K distribution o f ~7,200 stars with good photometry from the SDSS DR5 spectroscopic 
sample (jWest et al.ll2008l ). Using this sample, we derive the mean and standard deviation oi g — K 
for each spectral type, which allows us to estimate the spectral types of the ug-OYAy sample. 
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Fig. 3. — The distribution of M dwarfs sorted by the photometrically estimated spectral type for 
the total photometric sample (50, 130 stars), the total spectroscopic sample (4, 005 stars), the active 
stars (219 stars), and the inactive stars (3, 786 stars). 
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Fig. 4. — An example light curve with 46 observations in ugriz for an active M6 dwarf which exhibits 
both a large and small flare in the Stripe 82 dataset (highhghted with star symbols). The magnitude 
enhancements for the large flare are Amag= 5.50, 3.07, 2.34, 0.77, 0.29 in ugriz, respectively. Note 
that the small flare shows negligible enhancements in i and z. The inherent large variations in the 
u-band necessitate our requirement for a simultaneous enhancement in the gi-band. The panels to 
the right zoom in on 13 observations (over 78 days) around the small flare. The solid line indicates 
the quiescent magnitude whereas the dotted line shows the 0.7 magnitude flare detection threshold. 
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Fig. 5. — Example light curves in U and B of time-r esolved broadband photometry during a ~ 0.6 
U-hand peak magnitude flare (IKowalski et al.ll2009l . in prep). The time-sequence of the SDSS u 
(blue) and g (green) fllters are shown as the hashed regions in the unlikely scenario that the w-band 
integration caught the flare exactly at the peak. The time sequence introduces uncertainty in the 
interpretation of the SDSS flare emission, especially for small, short-lived flares. 
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Fig. 6. — The distribution of the flare variabihty index for our candidate flare distribution (red 
curve; > 0, Fuj — F^^g > 0, and Fgj — Fg^q > 0) compared to the distribution of < (blue 
curve) shown for the same number of epochs (~770,000). The former comprises both candidate 
flares and random night-to-night photometric variation, while the latter results purely from random 
photometric variation. The vertical line indicates the threshold value of = 100 which provides 
< 10% false discoveries, and gives 673 candidate flares (see text). 




Fig. 7. — An example of the red leak in SDSS w-band images for two red stars (both g — r = 1.44) 
during a non-flaring epoch (run = 6590, camcol = 2, field = 218) at a high airmass of 1.64. These 
represent the worst cases of red leak encountered, with nearly complete separation of the red lead 
and true stellar imag es. A line at the parallactic angle with length equal to the displacement of 
~2.3" predicted from iFilippenkd (119821 ) is overplotted in red next to each star. If all the red leak 
flux were added to the u flux and not deblended, the red leak would contribute approximately 0.5 
magnitudes to the photometry for each star. 
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Fig. 8. — Top - The flux enhancement in the u filter vs. for $„g > 0, F„ .,■ — Fu g > 0, and 
Fgj — Fg^g > 0. The vertical line corresponds to critical threshold = 100 and the horizontal 
line corresponds to a Au mag = 0.7. The 271 flaring epochs are shown as colored circles, while the 
observations which are cut out by flags or image inspection are left as black circles. The contours 
correspond to density levels of (5, 10, 25, 75, 150, 300, 600, 1200) in bins of A^^^ = 5, Ae„ = 0.1. 
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Fig. 9. — Top, left - The flaring fraction for the entire sample, binned according to the photometri- 
cally estimated spectral type. The numbers indicate the number of flares in each spectral type bin. 
Top, right - The flaring fraction for the ugriz-onlj sample using i — z color bins, showing that using 
either color or spectral type isolates the same groups of stars. Bottom, left - The flaring fraction 
for the active stars with SDSS spectra only, binned according to the spectral type determined from 
each spectrum. The flaring fractions for the active stars are significantly larger than for the entire 
sample. Bottom, right - The fraction of epochs that are large flares, using several values of hujiare / 
\jBoi to define a large flare. Although it is easier to see low-luminosity flares in the late- type dwarfs 
(contrast effect), the high-luminosity flares can be detected at all spectral types. The number of 
flares in the (early, mid, late) type bins are (33, 29, 46), (24, 14, 35), and (16, 8, 25) for the triangles, 
circles, and squares, respectively. 
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Fig. 10. — Top - The flaring fraction for different spectral type bins as a function of \Z\ distance 
from the plane. The MO-Ml stars are plotted in bins of 0-200 pc, 200 - 400 pc, and 400 - 900 pc; M 
-M3 stars are plotted in bins of - 100 pc, 100 - 200 pc, and 200 - 300 pc; M4-M6 stars are plotted 
in bins of - 50 pc, 50 - 100 pc, and 100 - 200 pc. Note the differences in volume sampled for the 
different spectral types (due to the u apparent magnitude limit). Also, note the sharp decrease in 
flaring fraction as older populations are sampled farther from the Galactic plane. Bottom - The 
cumulative distributions of \Z\ distance for all epochs (black squares), active epochs (red triangles), 
and flare epochs (blue asterisks), in bins of 50 pc. The most distant flare occurs on a star that is 
~850 pc below the plane, while the most distant star is located at ~1500 pc below the plane. 
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Fig. 11. — Top - The distribution of m- magnitude enhancements during flares. The dotted lines 
show the cumulative distributions. Later-type (lower mass) stars have the most detected flares and 
generate very large magnitude enhancements (due to their smaller quiescent luminosity). When 
converted to luminosities (bottom), the flares on the higher mass stars are the most luminous. 
Both of these results are consistent with previous flare observations. The errors in luminosities are 
~50%, and we show a typical error bar. 
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Fig. 12. — The flaring frequency for active M4-M5 stars from Stripe 82 c ompared to the a verage 
?7-band flare frequency of the very active, nearby dM4.5e star, YZ CMi (ILacy et al.lll976l ). The 
good agreement between the two distributions demonstrates that flaring frequencies obtained from 
observing a large number of stars at a low cadence is consistent with observing a single flare star 
continuously. The SDSS Stripe 82 distribution contains 18 flares. 
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Fig. 13. — The flaring fraction for 3 divisions of Galactic latitude (top) and longitude (bottom). 
The flaring fraction decreases for hnes of sight with larger Galactic latitude, probably because a 
smaller fraction of the sample is close to the Galactic plane (and therefore hkely to be younger 
and magnetically active). There is no obvious trend with Galactic longitude. Note that the flaring 
fraction is a normalized quantity, so number density variations have been removed. 
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Table 1. M Dwarf Colors and Spectral Types 



Spectral Subtype 



MO 


(0.55,0.32) 


Ml 


(0.78,0.43) 


M2 


(0.95,0.52) 


M3 


(1.15,0.62) 


M4 


(1.41,0.76) 


M5 


(1.87,1.03) 


M6 


(2.09,1.14) 


M7 


(2.49, 1.35) 



0.0159 


0, 


.0039 


0.0039 


0, 


.0083 


0.0135 


0, 


.0051 


0.0051 


0, 


.0055 


0.0163 


0, 


.0061 


0.0061 


0, 


.0058 


0.0204 


0, 


.0076 


0.0076 


0, 


.0076 


0.0213 


0, 


.0086 


0.0086 


0, 


.0085 


0.0183 


0, 


.0084 


0.0084 


0, 


.0062 


0.0194 


0, 


.0072 


0.0072 


0, 


.0054 


0.0270 


0, 


.0107 


0.0107 


0, 


.0095 



Note. — The best fit parameters at each M dwarf spec- 
tral type for each 2D Gaussian distribution of r — i, i — z 
colors , determined from 20,000 M dwarfs in the lWest et al.l 
ll2008h DR5 spectroscopic sample with well-measured col- 
ors and well-determined types. S is the covariance matrix 

and is of the form S = ( ''r-i,i-z \ r^^^ 

probability for each spectral type is evaluated as 'p{x) = 
2^||s||i/2 exp{-^{x - fi)'^T,-^{x - fj.)), M = ^^"^ j are 

the mean colors for a given type, and x are the r — i, i — z 
colors for a given star. 
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Table 2. Flare Selection Criteria 



Criterium 



ugriz 



ug 



Spectroscopic{SDSS) 



Justification 



Clean Sample 



^u, quiet ^ 
' ^g, quiet ^ 
; 100 
Au mag > 0.7 

SDSS Photometry Flag Cuts 
Image Check 



39,633 (1,987,253) 
39,629 (1,031,706) 
39,620 (612,728) 

321 (463) 
243 (296) 
201 (231) 
195 (225) 



10,497 (524,981) 
10,497 (272,249) 
10,484 (158,550) 

132 (210) 
60 (74) 

44 (49) 
41 (46) 



4,005 (197,867) 
4,004 (102,582) 
4,004 (61,107) 

36 (51) 
29 (38) 
... ( ... ) 

23 (29) 



M star colors, no WD-dM 
^ug < from noise in u and g 
Eliminates possible eclipses 

FDR 
Red Leak 



Note. — The criteria for selecting flares are shown, separately for the ugriz sample, the tig-only sample and the spectroscopic sample 
(those stars with SDSS spectra). The number of stars is given, and the number of epochs is shown in parentheses. The final flare sample 
contains 236 stars with 271 flaring epochs. 
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Table 3. APO 3.5m DIS (R300) Spectral Sample 



SDSS ID 


Spectral Type 


Ha EW (A) 


Notes 


Activity 


Dates 


J2Uoy4U.oo+UUioiO. ( 


Mi 


1 A 




y 


OoOozO 


Ton /1 1 Q Q o'7 nn i n cc o cr 
JzU41oo.o /-UUiyoo.o 


Mo 


A '7 
4.7 




y 


AOAOI n 


Ton/ii/i/i in 1 r\f\ A f\ A Q o 
J2U4i44.iU+UU4U4o.z 


TV /T /I A (f /I 

M4, M4 


O.'l, 1.5 




n, n 


Uoiiu/, uouyiy 


JzU445z.77+UiUo4D.y 


A /T 1 A TO 
Mi, MZ 


0.3, 0.9 




var 


A01 1 A'7 AQAAI A 

Ooii07, Oooyiy 


JzU472o.oU+UU4i^O.U 


AJO 

Mz 


3.3 




y 


AQAAI A 


Ton /I TOO 1 o 1 nnoon^? 

J 204728. IzH-OUzzOo.o 


M4 


8.9 




y 


AOAA1 A 


Toncmn oq mnoQn n 
J2UoU2U.oo-UiUooy.y 


ATQ 

Mo 


3.5 




y 


nonm n 

uouyiy 


Tonf;Qnn nc; nnuooi c 


Mo 


Q Q 
O.O 




y 


nonm n 
UoUyiy 


Tnncrr'oo cri i nm /i 1 o cr 

JzU5Doz.5i+UUi4i^i.o 


Ml 


2.2 




y 


nonm n 

uouyiy 




Mo 


3.8 




y 


nonm n 

Uouyiy 


Toinnoo Tn nncoon i 

j2lUUzo.7y-UU5ooU.l 


MO 


0.8 




w 


nonm n 

uouyiy 


JziUUo4.Uz-UU0ol4.0 


Mb 


6.3 




y 


nonm n 

uouyiy 


ToinU/iQ CO 1 nnnono i 


A /TQ 

Mo 


7.0 




y 


nonm n 

Uouyiy 


T01 1 1 1 n OCT 1 nncrno'? i 


M4 


A C 

4.0 




y 


nonm n 

Uouyiy 


701 1 1 QQ f\A 1 nnOI Q'7 o 


M4 


0.9 




y 


nonm n 

Uouyiy 


T01 1 Q1 1 etc nnncnQ o 


M4 


4.1 




y 


nonm n 

Uouyiy 


Jzil414.Uo-UUuOo4. / 


Air; 
Mo 


5.0 




y 


nonm n 

uouyiy 


7m 1 /( cT/i CO nnnoi c q 
J z i 14o4. OZ-UUUz lo . o 


M4 


5.3 




y 


no 1 1 n'7 
UollU7 


T01 1 oo'7 nn 1 nnnO/( q c 


A/TO A4'n A^"! 

MZ, MU, Mi 


1.6, 1.4, 0.8 




y, y, w 


noiin'7 nonooc noioi/i 
UollU7, UoUozo, Uolzl4 


T01 01 QQ on ni 1 1 Q'7 A 

Jzizioo.oy-Uiiio / .4 


M4 


5.5 




y 


nonm n 

uouyiy 


jziz / 4y. /y+uuzDU / .z 


Mo 






y 


UoiiU ( 


JzloliZ.Uo-UUzlo4.U 


M4 


4.0 




y 


no 1 1 n'7 
UollU7 


JzlozUo.z / +UU4UzD.y 


Mo 


7.0 




y 


nonm n 

Uouyiy 


J2iooi7.y4^UU57zo.i 


ATI 

Mi 


1.6 




y 


nonm n 

Uouyiy 


T01 dnon 1 nnQQon a 
Jzl4Uzy.z4-|-UUoo2U.O 


Mo 


1.0 




w 


nonm n 

uouyiy 


701,110/1 CI nnooon c 


"KX A 

M4 


2.0 




y 


no 1 oi 
Uoizi4 


Jzl4zzb.iO+UiUo4Z.o 


Mi 


0.5 




w 


nonm n 

Uouyiy 


T01 /lO/ic; cT'T nncT/iQi^ v 
JZi4Z4o.O /-UUo4oD. ( 


1\ /Tl 

Mi 


1.2 




y 


nonm n 

Uouyiy 


T01 /ivcci tin nnQOQV a 
Jzl4/ol.Dy-UUoZo 1 A 


M4 


3.9 




y 


nonm n 

Uouyiy 


TmcT/iQn oo 1 nno i i T n 

Jzl54oU.oo+UUzll7.U 


M4 


3.3 




y 


nonm n 

uouyiy 


Toonc; /( Q on nnnQ co c 


ATI A(T1 
Mi, Mi 


0.0, 0.2 




n, n 


nonooc nonm n 


Toonnon on nncr'voo o 

JizUyiU.iy-UUoYjs.o 


M2 


2.4 




y 


no 1 01 y( 

Uolzl4 


TOOT Tor' "1 n 1 nnoncro o 

Jzzl7ZD.lU+UUoUoZ.z 


M3 


2.9 




y 


no 1 m o 

UolUlo 


Tooom o on nnQ/ion n 

jzZzUiz.oy-uuo4zy.u 


A /Tl 
Mi 


1.3 




y 


nonoo c: 
UoUozo 


Tooocico cc/i 1 nnno i t; 


M4 


A 1 

4.i 




y 


nonm n 

uouyiy 


TOOQ'71 1 ni 1 nncQnfi c 


ATO 

Mz 


2.6 




y 


nonm n 

uouyiy 


JZZo < ol.OO-t-UUo4Zo.O 


IVlO 


o.y 




y 


UoiUio 


TOO crz?y(i oT 1 minrrn n 

J225o41.o7+UllUoU.y 


M5 


2.8 




y 


nonoo c 

080825 


J230001.05+010454.1 


M5 


5.6 




y 


081107 


J230314.24-003110.7 


M2 


2.7 




y 


080825 


J231109. 45+002515.4 


M3 


4.2 




y 


081107 


J231810.79-004512.3 


M2 


3.4 




y 


080825 


J232516.19-001600.9 


M2 


1.8 




y 


080825 


J233810.08-004450.3 


M4 


4.3 




y 


080825 


J233952.93-001005.3 


M5 


5.8 




y 


080825 


J234612. 94+002507.3 


Ml 


0.0 




n 


081107 


J004145.91-002832.7 


M5 


6.0 




y 


080825 


J005847.77+011015.5 


M4 


4.5 




y 


080825 


J010052.58+010233.1 


MO 


-0.2 




n 


081107 
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Table 3 — Continued 



snss TD 




Fry F,W ("Al 




III V 1 Lr y 




J011745. 28+005929.4 


Ml 


1.8 




y 


081013 


J012531. 29-005359.0 


M3 


2.7 


_ 


y 


080825 


J013851.46-001621.7 


M6 


12.0 


_ 


y 


081013 


J020056. 76+004149.1 


M4 


9.2 


_ 


y 


081013 


J020327.92-000442.1 


M3 


3.1 


_ 


y 


081013 


J022112. 48+003225.7 


M4,M3 


0.7,0.8 


_ 


n,n 


080825, 080919 


J022214.35-003217.5 


M2, M2 


0.7,0.6 




w, w 


081013, 081124 


J023452.37-010031.4 


M3 


3.3 


- 


y 


081107 


J023538.57-000422.2 


M4 


6.0 




y 


081107 


J023547.98-005822.4 


M3 


2.1 




y 


081107 


J023550.72+010136.8 


M3 


2.3 




y 


081124 


J024554. 14-010445.1 


M3 


6.0 




y 


081124 


J025747.56+004731.2 


M4 


3.8 




y 


081124 


J032605.67-002402.6 


MO 


2.0 




y 


081124 


J033017.55+005359.4 


K7 


-0.5 




n 


081124 


J033734.14-011104.5 


M2, M3 


-0.6, 0.0 




n, n 


081013, 081124 


J033735. 25+005859.7 


M2 


2.4 




y 


081124 


J205059. 98+005136.4 


M4 


2.4 


ug-onlj 


y 


080919 


J205257.83-000335.9 


M3 


3.3 


ug-only 


y 


080919 


J215455. 18+011414.5 


MO 


6.0 


ug-only 


y 


081214 


J215517.42-004547.9 


M5 


11.6 


ug-only 


y 


081214 


J221712.10-000815.4 


MO 


-0.4 


ug-only 


n 


081214 


J001208. 37+002559.1 


M4 


5.7 


ug-only 


y 


081214 


J005903.60-011331.1 


M2 


2.6 


ug-only 


y 


081214 


J011307.16-002326.1 


M3 


3.4 


ug-only 


y 


081214 


J012904.68+003058.4 


M5 


2.3 


ug-only 


y 


081214 


J020911. 90+004432.4 


Ml 


1.2 


ug-only 


y 


081214 



Note. — y - active; w - weakly active; n - inactive; var - variable. Activity status was determined using the criteria in I West et al.l 
l|2004l .l 20081 ) and types were determined automatically with the Hammer facility (Covey et al. 2007). The automatically assigned activity 
status and type were adjusted by eye as needed, and EW's were recalcula ted using accounting for the r adial velocity shift. Errors in the 
EW are are typically ~0.2 A, and were calculated using Equation AlO in IChalabaev fc MaillardI l ll983h . 
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Table 4. Spatial Flare Rates 



Spatial Division 


Mean # flares hr ^ sq deg ^ 


Flaring Fraction 


23° < |fe| < 45° 


2.0 


1.38x10""' 


45° < |6| < 55° 


0.9 


9.9x10-5 


55° < \b\ < 64° 


0.7 


8.0x10-5 


45° < i < 90° 


1.5 


1.11x10-* 


90° <l< 135° 


0.7 


8.3x10-5 


135° < I < 191° 


0.9 


1.16x10-* 



Note. — We present the raw spatial rate (mean number of flares hr— ^ 
sq dcg-^) and the flaring fraction for different ranges of Galactic latitude 
and longitude. When calculating the raw spatial rate, we only consider 
the 1 X 1 sq deg bins with at least 50% SDSS coverage. 



-40- 



Table 5. Flare Rate Model Parameters 



Spectral Subtype Bin 




f active 


fa, flare 


/inactive 




limiting log Ljy 


Predicted # flares* 


Observed # flares 


MO-Ml 


35,000 


0.03 


1.0 


0.97 


0.02 


29.5 


20-66 


38 


M2-M3 


13,000 


0.06 


1.0 


0.94 


0.1 


28.7 


68-150 


83 


M4-M6 


2,100 


0.30 


1.0 


0.70 


0.4 


28.0 


193-233 


150 



Note. — We present the measured and derived values for the flare rate model given by Equation 4 in §4.6.1. 

*We adopt 1% and 10% for fijiarei the fraction of inactive stars that are capable of flaring, which give a range of values for the 
predicted number of flares. 
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